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We deduce the effects of quantum interference on the con- 
ductance of chaotic cavities by using a statistical ansatz for 
the S matrix. Assuming that the circular ensembles describe 
the S matrix of a chaotic cavity, we find that the conductance 
fluctuation and weak-localization magnitudes are universal: 
they are independent of the size and shape of the cavity if the 
number of incoming modes, N, is large. The limit of small 
N is more relevant experimentally; here we calculate the full 
distribution of the conductance and find striking differences 
as A^ changes or a magnetic field is applied. 

The effect of quantum interference on transport 
through microstructures has been intensively investi- 
gated and is one of the main subjects of mesoscopic 
physics B. For diffusive transport in disordered struc- 
tures, both microscopic perturbative and macroscopic 
random matrix theories give a good account of the phe- 
nomena. In the latter case |2|, one assumes that the 
transfer matrix for an ensemble of such disordered mi- 
crostructures can be chosen from a simple statistical en- 
semble with only symmetry constraints applied. The suc- 
cess of this theory is perhaps the best theoretical demon- 
stration of the ubiquity of mesoscopic interference effects. 

More recently, interest has focused on quantum trans- 
port in ballistic microstructures — structures in which 
impurity scattering can be neglected so that only scat- 
tering from the boundaries of the conducting region is im- 
portant 0]. Quantum interference effects in such struc- 
tures depend on the nature of the classical dynamics 
|S-@|, in particular whether it is regular or chaotic ||. 
Recent experiments have studied transport in such bal- 
listic structures Jl0|-^3[ and have detected a difference 
between nominally regular and chaotic structures |lC[ ]. 

The theoretical work on this subject |j]-[7| has con- 
centrated on cither numerical quantum calculations or 
semiclassical theory. On the other hand, it has been pro- 
posed |8| that chaotic scattering in the quantum regime 
JLlj should be described by a random matrix theory for 
the S-matrix. The emphasis in both that work and very 
recent work on the iS-matrix of disordered structures 
is on the properties of the eigenphases of S. The eigen- 
phases, however, are not directly connected to the trans- 
port properties because they involve both reflection and 
transmission. In contrast, in this paper we derive the 
implications of such a random 5-matrix theory for the 
quantum transport properties and provide numerical ev- 
idence that this theory applies to the class of ballistic 



microstructures investigated experimentally. In this way 
we obtain experimentally accessible predictions for the 
quantum transport properties of chaotic billiards. 

A quantum scattering problem is described by its S- 
matrix. For scattering involving two leads (see Fig. |l|) 
each with N channels and width W, we have 

»-(;$) w 

where r, t are the N x N reflection and transmission ma- 
trices for particles from the left and r' , t' for those from 
the right. In terms of S, the conductance is JTg[ 

G = (e 2 /h)T = (e 2 /h)Tr{t^}. (2) 

Because of current conservation S is unitary, SS* = I, 
and in the absence of a magnetic field it is symmetric 
because of time-reversal symmetry, S = S T . 

We concentrate on situations where the statistics of the 
scattering can be described by an ensemble of S-matrices 
that assign to S an "equal a priori distribution" once the 
symmetry restrictions have been imposed. In particular, 
the possibility of "direct" processes caused, for instance, 
by short trajectories and giving rise to a nonvanishing 
averaged ^-matrix (lj]||, is ruled out. The appropriate 
ensembles are well-known in classical random matrix the- 
ory p8| and are called the Circular Orthogonal Ensemble 
(COE, f3 = 1) in the presence of time-reversal symmetry 
and the Circular Unitary Ensemble (CUE, j3 = 2) in the 
absence of such symmetry. These ensembles are defined 
through their invariant measure: the measure on the ma- 
trix space which is invariant under the appropriate sym- 
metry operations. To be precise |Q, d^i(S) = dfj,(S'), 
where S' = UqSVq, and Uq, Vq are arbitrary fixed uni- 
tary matrices in the case of the CUE with the restriction 
Vq = Uq in the COE. Numerical evidence for the valid- 
ity of this random matrix theory for describing quantum- 
chaotic scattering can be found in Ref. || . 

Perhaps the most widely studied mesoscopic trans- 
port effects are the magnitude of the conductance 
fluctuations — how much the conductance varies as a 
magnetic field or gate voltage is applied — and the size of 
the weak-localization correction to the average conduc- 
tance at B = M. We therefore start by deriving (T) 
and var(T) where we use an integration over the invari- 
ant measure as the average. Such integrals have been 
evaluated previously , and we find that 
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Performing the trace over channels in Eq. (|^) , we obtain 
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where the limit is as N — ► oo. 

We make several comments concerning these results. 
(1) Previously, semiclassical theory and numerical calcu- 
lations suggested that the weak-localization correction, 
(T) — N/2, and the magnitude of the conductance fluc- 
tuations, var(T), are independent of the size of the sys- 
tem for chaotic billiards ||] . This is the analogue of the 
"universality" of the conductance fluctuations in the dif- 
fusive regime jjj . Since the number of modes is propor- 
tional to the size of the system (N — int[kW / -k]) , our 
N — > oo results show that the conductance fluctuations 
and weak-localization are universal within the random S- 
matrix theory. (2) In the large N limit, var(T) in the 
presence of time-reversal symmetry is twice as large as 
in the absence of symmetry, as for diffusive conductance 
fluctuations. This demonstrates the universal effect of 
symmetry. (3) Both quantities show some variation in 
the small N regime typical of the experiments [p~o| — P~2]| - 
Note, for instance, that the ratio of var(T) in the pres- 
ence and absence of symmetry is not 2 for small TV. (4) 
The values obtained in the N — > oo limit are the same 
as those obtained from a random matrix theory for the 
Hamiltonian EQ] in which one assumes that the Hamilto- 
nian of the billiard is described by the Gaussian Ensem- 
bles and finds the conductance by coupling the billiard to 
leads in a random way. Quantitative agreement between 
S'-matrix and Hamiltonian random matrix theories has 
been noted in the past (I?]] but is not fully understood. 

The predictions of the random matrix theory are com- 
pared to the conductance of a stadium billiard in Fig. [j] 
computed using the method of Ref. [pi"[ . I n these cal- 
culations, the S-matrix varies as a function of energy 
because of the resonances occurring in the cavity. We 
estimate that the resonances are moderately overlapping 
for N = 1 and that the width to spacing ratio increases 
linearly with N. The basic assumption (ergodic hypothe- 
sis) is that through these fluctuations S covers the matrix 
space with uniform probability. This should apply to bil- 
liards in which the effect of short non-chaotic paths is 



minimized. We therefore use a stadium billiard in which 
(1) a stopper blocks any direct transmission between the 
leads, (2) a stopper blocks the whispering gallery tra- 
jectories which hug the half-circle part of the stadium, 
and (3) the stadium is asymmetrized to break all reflec- 
tion symmetries. We obtain excellent agreement between 
the energy averages found numerically and the invariant- 
measure ensemble averages introduced above. In particu- 
lar, both the variation at small N and the ratio of var(T) 
in the presence and absence of time-reversal symmetry 
are verified in the billiards. 

Motivated by this good agreement, we consider more 
detailed predictions of the random S'-matrix theory: we 
derive the full distribution of T for small N and the statis- 
tics of the eigenvalues of W, denoted {t}. We obtain 
these results by expressing the invariant measure in terms 
of a set of variables that includes the {r}. Any unitary 
matrix of the form in Eq. (Q) can be written as 0] 
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where r stands for a N x N diagonal matrix of the eigen- 
values {t} and the v^' are arbitrary unitary matrices ex- 
cept that = (v^) T and = (v^) T in the presence 
of time-reversal symmetry. Now it is a general property 
of measures on vector spaces |^3| that a differential arc- 
length written in the form da 2 — ^2 ab g a bdx a dx b implies 
that the volume measure is d[i(V) = y/det{g) ]J a dx a . 
In our case the differential arc-length is simply da 2 = 
Tr{dS"<dS}. Substituting for S the form in Eq. (@), one 
finds (fi = 1,2) 
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where the joint probability distribution of the {r} is 



P2({r})=C 2 Y[\r a 
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dfj,(v^) denotes the invariant, or Haar's, measure on the 
unitary group [24|, and Cp are N dependent normaliza- 
tion constants 25]. 

The distribution of T = Y^ a =i Ta f° uows from Eq. 
(||) by integration over the joint probability distribu- 
tion. This can be carried out for small N; for instance, 
in the trivial case N = 1, w(T) = 1 for the CUE and 
w(T) = 1/(2VT) for the COE. For N = 1 - 3 the w{T) 
derived from the random matrix theory are plotted in 
Fig. ^| and compared to numerical data for billiards. 
Note the dramatic difference between the CUE and COE 
w(T) in the single mode case, and the difference within 
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each ensemble between the N = 1 and N = 2 cases. The 
results for TV = 3 are close to a Gaussian distribution 
denned by the two moments given in Eqs. (^). 

The agreement between the numerical and theoreti- 
cal results in Fig. is very good in terms of both the 
difference between COE and CUE and the dependence 
on N ^6|. These effects should be observable in exper- 
iments, which are typically done in the few mode limit, 
and would provide a clear test of the applicability of ran- 
dom S- matrix theory to experimental microstructures. 

Though not experimentally accessible, the {r} are the- 
oretically interesting because of their fundamental rela- 
tion to the conductance. We obtain further information 
by writing the joint probability density in the form 

Pp({r}) = Cp exp{-/3[^ In \r a - n\ + £ Vp{r c )]} (7) 
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with 14(r) = and V\{t) = |lnr. This is exactly the 
form of the joint density in the global-maximum-entropy 
approach to quantum transport in disordered systems g] 
and in the Gaussian ensembles p8|| ■ Many statistics of 
such distributions are known asymptotically as N — > oo 
jl8],^7 - 30 1 . For instance, the known form of the asymp- 
totic two-point correlation function [^7|,[50) can be used 
to obtain the asymptotic value of variT): 
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This result agrees with Eq. (^g) and is independent of 
the potential V(r). Therefore the asymptotic value of 
var(T) is the same for a large class of random matrix 
theories, another aspect of the "universality" of conduc- 
tance fluctuations |^9) . 

In the CUE case, the joint density in Eq. (Jfy is of a spe- 
cial form suitable for the random matrix theory method 
of orthogonal polynomials fl^,||]. In this case, V = 
and the eigenvalues are restricted to [0, 1] so the Legen- 
dre polynomials are appropriate fjl|| . In terms of these 
polynomials, the exact eigenvalue density and two-point 
correlation function are piLp7| 
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p 2 (r,r') = P (t)S(t-t')/N (9) 
~[P N {a)P N ^{a') ~ P N ^{a)P N {a')] 2 /\Q{t - t') 2 

where a = It — 1. Using the asymptotic expansion of the 
Legendre polynomials as N — > oo and some smoothing, 
one finds that p(r) — > N/tt^t(1 — r) and recovers 
the expression in Refs. [^j and |30| for the asymptotic 
two-point correlation function. Previous work has shown 
that the statistics of the eigenvalues {t} follows that of 
the Gaussian unitary ensemble in the large N limit [BlJ . 



Comparison of the predictions of the last two paragraphs 
with numerical results further supports the applicability 
of the COE and CUE to ballistic cavities. 

In summary, we have derived the consequences for 
quantum transport of the assumption that the S'-matrix 
of a chaotic cavity follows the circular ensembles. We 
have shown that the magnitude of both the conductance 
fluctuations and the weak-localization in a chaotic mi- 
crostructure are universal in the large N limit. The small 
N limit is most relevant experimentally, and here we find 
deviations from the asymptotic behavior (Fig. |l|) as well 
as a striking dependence of the full distribution of T on 
both N and magnetic field (Fig. ^). In closing, we em- 
phasize that we have neglected the "direct" scattering 
due to short paths ((5) = 0); since such scattering is 
important in many chaotic cavities, the effect of these 
processes on quantum transport remains an important 
open question, which, in principle, could be investigated 
using the information-theoretic model of Ref . (l7) . 

We thank R. A. Jalabert for several valuable discus- 
sions. One of us (PAM) appreciates the kind hospitality 
of the Wissenschaftkolleg zu Berlin where part of this 
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Note added — While preparing this work for publica- 
tion we received a preprint by C. W. J. Beenakker, R. 
A. Jalabert, and J.-L. Pichard which contains some over- 
lapping material. We thank J.-L. Pichard for sending us 
this preprint. 
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FIG. 1. The magnitude of the weak-localization correc- 
tion [panel (a)] and the conductance fluctuations [panel (b)] 
as a function of the number of modes in the leads, N. The nu- 
merical results for B — (squares with statistical error bars) 
agree with the prediction of the COE (diamonds connected by 
dotted line), while those for B / (triangles) agree with the 
prediction of the CUE (diamonds connected by dashed line). 
The inset shows a typical ballistic cavity used. The numer- 
ical results involve averaging over (1) energy at fixed N (50 
points), (2) 6 different cavities obtained by changing the two 
stoppers, and (3) 2 magnetic fields for (BA/<j>o = 2, 4 

where A is the area of the cavity). 

FIG. 2. The distribution of the transmission intensity at 
fixed N = 1, 2, or 3 both in the absence (first column) and 
presence (second column) of a magnetic field. The numerical 
results (plusses with statistical error bars) are in good agree- 
ment with the predictions of the circular ensembles (dashed 
lines). Note in particular the striking difference between the 
N = 1 (first row) and N = 2 (second row) results and be- 
tween the B — and B / results for N = 1. For N = 3 the 
distribution approaches a Gaussian (dotted lines). The cavi- 
ties used are the same as those in Fig. 1; for B ^ 0, BA/<j>o = 
2, 3, 4, and 5 were used. 
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